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ABSTRACT 

We present a mid-infrared spectroscopic data cube of the central part of 30 Doradus, observed 
with Spitzer^s IRS and MIPS/SED mode. Aromatic dust emission features and emission hnes 
from molecular and atomic hydrogen are detected but not particularly strong. The dominant 
spectral features are emission lines from moderately ionized species of argon, neon, and sulphur, 
which are used to determine the physical conditions in the ionized gas. The ionized gas excitation 
shows strong variations on parsec scales, some of which can plausibly be associated with individual 
hot stars. We fit the ionic line strengths with photoionization and shock models, and find that 
photoionization dominates in the region. The ionization parameter U traces the rim of the central 
bubble, as well as highlighting isolated sources of ionization, and at least one quiescent clump. 
The hardness of the ionizing radiation field Trad reveals several "hot spots" that are either the 
result of individual very hot stars or trace the propagation of the diffuse ionizing field through 
the surrounding neutral cloud. Consistent with other measurements of giant H II regions, log{U) 
ranges between -3 and -0.75, and Trad between 30000 and 85000K. 

Subject headings: HII regions - ISM: individual (30 Doradus) - infrared: ISM - Magellanic Clouds 
1. Introduction 
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The 30 Doradus region of the Large Magellanic 
Cloud (LMC) is an ideal laboratory in which to 
study the effect of massive star formation and its 
feedback on the circumcluster interstellar medium 
(ISM). Star formation processes in the Magel- 
lanic Clouds are potentially a template for the 
early universe, where small irregular galaxies were 
commonplace and the overall metallicity lower. 
In particular, the star formation rate measured 
relative to molecular mass may be high in the 
LMC compared to our Galaxy, but that measure- 
ment is complicated by interpretation of CO data 
in molecular clouds that are known to be more 
porous and deeply penetrated by ultraviolet radi- 
atio n, and an unc ertain C O-H2 conversion fa. ctor 
(e.g. IPoglitsch et al..,1995t [Bernard et all2008( ). 
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30 Doradus itself, the brightest star for- 
mation region in the LMC, contains several 
lO^ MfT) of molecular hydro gen traced by ^^C0(1- 



0) ([Johansson et alj 119981 ) , much of which is in 



two elongated clouds that form an arc or ridge in 
the center of the nebula (see Figure [1] for a multi- 
wavelength view and Figure [2] for the distribution 
of CO). That molecular material is probably only 
the remnant of the cloud that formed the cen- 
tral cluster and thousands of OB stars. Observa- 
tions of higher energy CO transitions suggest that 
this remna nt ridge of m olecular gas is quite warm 
and dense ( KimI l2007l l. and near- infrared obser- 



vations ha ve revealed that it i s actively forming 
new stars (iHyland et al.lll992t iRubio et aLlllQQSl: 
Maercker fc BurtonlbOOSh 

The dominant star-forming cluster of the 
30 Doradus H II region is NGC 2070, whose dense 
center, R136, is usually considered the nearest su- 
per star cluster (SSC). R136 contains 39 03 stars, 
a total stellar mass of ~5 x IC^Mq within 2.5 pc, 
and stellar densities exc e eding 5 x IO^Mqpc""^ 



(|Hunter I Il999l : IWalborn I Il99ll ). Extragalactic 



SSCs are the hosts of one of the most extreme 
modes of star formation in the universe, may de- 
velop into globular clusters, andplay a key role in 
galaxy formation ( Johnson 20041 ) . At a distance 
of 50±2.5 kpc (see discussion of th e LMC distance 
and uncertainty in ISchaefer 2008 ). 30 Doradus is 
close enough to study at the parsec-scale resolu- 
tion that is required to understand the formation 
and feedback effects of individual stars, and a 
prime target for detailed study of the same mech- 
anisms which operate in more distant and massive 
SSCs. 

Observations at mid-infrared (MIR) wave- 
lengths offer several advantages for studying star 
formation and its interaction with circumclus- 
ter dust and gas. MIR (here roughly defined as 
5-50/im) continuum emission is generally domi- 
nated by radiation from very small dust grains 
(VSGs); the broad shape of the continuum is 
sensitive to the VSG size and temperature dis- 
tribution, and hence indirectly to the radiation 
intensity in photon-dominated regions, and to the 
destruction of smaller grains as would be expected 
above 2GGGK and in HII regions. Broad (~l/im) 
dust emission features are also present between 3 
and 19/im. These previously named "unidentified 
infrared bands" result from distortion (bending, 



stretching) modes of aromatic molecules contain- 
ing tens to hundreds of carbon atoms. (Most 
features have been attributed to polycyclic aro- 
matic hydrocarbons, PAHs, but unique astrophys- 
ical identification is a work in progress). Analy- 
sis of relative strengths of PAH features reveals 
their size and ionization state which is expected 
to change in intense radiation fields. Extinction 
by dust in the infrared is low compared to other 
wavelength regimes. Observations in the MIR 
can pierce cold molecular clouds and reveal the 
star-forming regions that they shroud. There are 
two major bands of absorption by silicate dust 
at 9.7 and 18/im, whose shape and strength pro- 
vide further diagnostics of the dust in the region. 
Mid-infrared (MIR) spectroscopy o f the entire 
30 Doradus region with ISO-SWS (jSturm et al ' 



200n[ ) revealed a continuum-dominated spectrum 



with very weak silicate absorption and also only 
modestly strong aromatic emission features. The 
4 spectra taken in 30 Doradus with ISOPHOT 
( Vermeil et al. 20021 ) similarly show low ratios of 
PAH strength to IR continuum. 

The most recent and sensitive MIR contin- 
uum observations of the Magellani c Clouds were 



obtained as part of the Spitzer (jWerner et al 
2004) Legacy program "Surveying the Agents 
of a Ga laxy's Evolution" (SAGE), using IRAC 



f3-8um:lFazio et al-l liool and MIPS (24-160/im; 



Rieke et al.ll2004[ ). The goals of SAGE are to con- 
duct a detailed study of the dust processes in the 
ISM, to complete a census of newly formed stars 
in order to find the star formation rate in the 
LMC, and to study evolved stars with high mass 
loss rates in order to determine the rate of mass 
injection into the ISM (.M ei xner et al. 2006). Fig- 
ure [T] places the MIR emission from 30 Doradus in 
its multiwavelength context. 8/im emission traces 
the same arc-shaped ridge seen in the optical, de- 
lineating the edge of a bubble probably blown by 
R136 and filled with hot X-ray emitting plasma. 
The remaining molecular material in the region is 
also located in that ridge. 

This paper presents a new spatially filled spec- 
tral cube of 30 Doradus with the low-resol ution 
modules of Spitzer/IRS (jHouck et all l2004l) an d 



with theSEDmode of Spitzer/MIPS ( Rieke et al 



2004 iLu et all 12008) . The sensitivity of this 



dataset gre atly exceeds the pr evious ISOCAM- 
CVF data ([Madden et al.ll2006h . which only pro- 
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vided maps of the few strongest ionic lines, but 
tentatively showed a gradient of decreasing exci- 
tation in the [Ne III]/[Ne II] ratio as a function of 
distance from R136, and signs of PAH destruction 
in the central regions of the nebula. 

The data are described in sections fj2l In par- 
ticular, H2.1\ and Appendix |A] discusses reduction 
and artifacts, §2.3l the line and feature fitting pro- 
cedure, and Appendix |B] quality assurance tests 
including a comparison with high-r esolution IRS 
data from iLebouteiller et aL 1 (|2007l) . In ^ we 
describe the general results - the spatial distribu- 
tion of emission lines and their ratios. §4]describes 
the derived distribution of matter - electron den- 
sity and dust, and [J5]the distribution of radiation 
evident in the excitation of the gas. We compare 
the data with photoionization and shock models, 
and summarize and discuss implications in ij6l 

2. Observations 

The Spitzer^s Infrared Spectrograph has >100 
times greater spectroscopic sensitivity than the 
previous premier infrared observatory, ISO-SWS, 
and >10 times higher spatial resolution than the 
SWS aperture size. We used the four low resolu- 
tion modules of the I RS. The spectral resolution 
range s from 60 to 120 (jSpitzer Observer's Manual 7.1 
200# 1. with reliable wavelength coverage as noted 
in section IA.0.31 The four modules are the short 
wavelength / low resolution second order (SL2), 
short wavelength / low resolution first order (SLl), 
long wavelength / low resolution second order 
(LL2), and long wavelength / low resolution first 
order (LLl). Our data is divided into eighteen As- 
tronomical Observing Requests (AORs), each con- 
taining many BCD (basic calibrated data) frames, 
which are pipeline-processed images of the IRS 
chip. 

We observed 30 Doradus over five days in 
September of 2006, using Spitzer's spectral map- 
ping mode in order to obtain detailed spatial in- 
formation. The total amount of time on target 
was 74 hours. We used 3440 slit pointings, cover- 
ing 40.5 square arcminutes in the short- wavelength 
modules and 69.1 square arcminutes in the long- 
wavelength modules (see Figure [3]). Each slit over- 
laps half of the preceding one, and each row of 



slits overlaps half of the preceding row, so every 
point on the map was observed at four separate 
slit pointings. There were three repetitions per 
pointing for the SL observations, at 14 seconds 
each (for a total time of 168 s for each point on 
the SL maps), and four repetitions per pointing 
for the LL observations, at 6 seconds each (for a 
total time of 96 s for each point on the LL maps). 
The result is four large and dense spectral cubes, 
one for each module. 

We also made separate observations of a nearby 
region of comparatively blank sky, bracketing our 
target observations in time, to comprise back- 
ground observations and to better characterize 
rogue pixel response. The total time on back- 
ground was 44 minutes. 

2.1. Reduction of IRS data cube 

We used the basic calibrated data (BCD) 
from the Spitzer Science Center pipeline version 
S14.4.0. The main tool for the assembly and re- 



duction of our data cube was CUBISM ( Smith 



2007br ). which is designed for spatially dense IRS 
maps. The functions of CUBISM include tools to 
reduce the data (subtracting background, apply- 
ing a slit loss correction function, trimming the slit 
in the spatial direction, masking pixels that are 
fiagged by the pipeline, and applying algorithms 
to identify rogue pixels) and to extract spectra 
and maps from the data cube. 

Appendix [X] describes in detail our procedure 
for calibrating and adjusting the fiux levels such 
that the different spectrograph modules could be 
combined into smooth spectra. In brief, we used 
contemporaneous off-source background observa- 
tions to subtract the thermal background. We 
used measurements of the slit-loss correction func- 
tion provided by the Spitzer Science Center to 
undo that part of their point-source-based calibra- 
tion (since our observations are more similar to a 
uniform diffuse source than a point source). Fi- 
nally, we adjusted spectra from different modules 
to match each other. 

Appendix |B] describes our comparisons of high 
and low-resolution spectra in 30 Doradus to test 
the fidelity of our fitting process. We find that 
the strengths of even somewhat blended lines are 
recovered well in the low-resolution spectra. 



1 Available at; |http: / /ssc. spitzer. caltech.edu/documents/SOM/ 1 
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Fig. 2. — IRAC 8^m image (log scale). White 
contours: high -resolution 3cm continuum from 
Lazendic et al. 1 dioos). Contour levels: 6, 9, 12, 
15, 18, and 21 x 10~^ Jy/bm. Sources of inter- 
est marked here are listed in Tabled described in 
§3.3| and their spectra shown in Figure flTl 



Fig. 1.— The 30 Doradus nebula. Red: IRAC 
8/xm image (SAGE). Green: ESO B-band im- 
age. Blue: broad-band soft X-ra y image, 0.5 — 2 



keV (private correspondence and iTownslev et al 
20061) . All images are on a line ar scale. White 



contours: CO (1-0) emission ([Johansson et al 
[logs') . Black contour: a single level of 3cm radio 
emission, to guide the eye (see Figure [2|). Magenta 
mark: the star cluster R136, core of NGC 2070, at 
5''38"42^ -69°06'00". Cyan mark: the star cluster 
Hodge 301, at 5''38™16^ -69°04'00". Yellow mark: 
an IR p oint source in the molec ular cloud [JGB98] 
Dor-06 ([Johansso n et al.l[l998h . at 5''38'"31.63^ - 
69°02'14.6". This source is also marked in white 
in Figure [4| and discussed in 
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Fig. 3. — IRAC 8/xm image with the scope of cov- 
erage in each of the IRS modules shown. Red: 
LLl. Yellow: LL2. Blue: SLl. Magenta: SL2. 
Green: high resolution GTO apertures (see Ap- 
pendix [BJ. 
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2.2. Artifacts in IRS data cube 

There are two major artifacts which appear in 
a 2-D map of our data cube at any given wave- 
length. In all four modules, there are numerous 
faint stripes, one pixel wide, which cross the map 
in the direction that the slit scans. In the two 
SL modules, there is a broad intensity discontinu- 
ity and associated bright region in one area of the 
map. The affected region varies with wavelength. 
The faint stripes appear to be mainly caused by 
hot pixels on the IRS chip. As the slit was stepped 
across 30 Doradus, each hot pixel was "dragged" 
across the map, creating a bright stripe. They can 
be eliminated via two cleaning methods: wavsamp 
trimming and rogue pixel masking. The intensity 
discontinuity is caused by saturation in the peak- 
up image, and must be corrected by fitting a cor- 
rection factor to the regions of the chip between 
the SL modules. 

The Spitzer Science Center has defined a poly- 
gon, called the wavsamp, defining the active area 
to extract from the IRS chip for each module. In 
the case of 30 Doradus, the default wavsamp tends 
to be too generous in the spatial direction, includ- 
ing a few pixels where the spectral response is re- 
duced. By trimming it, we were able to eliminate 
some of the faint stripe artifacts. 

It is possible to determine which pixels on 
the IRS chip contribute to a given point in the 
data cube, using a CUBISM tool called backtrack- 
ing. Backtracking from the artifacts in the maps 
demonstrated that most, if not all, of the remain- 
ing faint stripes are caused by individual rogue 
pixels, whether hot or cold. In addition to mask- 
ing the pixels automatically flagged by the Spitzer 
Science Center pipeline, we employed CUBISM's 
automatic rogue pixel masking algorithms at both 
the global and record levels. At the global level, 
we masked any pixel which deviated by at least 
2.5 sigma from the median pixel level in at least 
35% of the records. At the record level, it was 
necessary to be much more conservative to avoid 
masking out real spectral features. We masked 
any pixel which deviated by at least 7.5 sigma in 
70% of its occurrences in the cube. The procedure 
missed some pixels which were clearly hot or cold, 
so we also manually marked a set of global and 
record level rogue pixels. In the resulting cubes, 
the faint stripes were greatly reduced. 



The other major artifact, the intensity discon- 
tinuity, is caused by saturated sources on the 
peak- up image (PUI), primarily by bright point 
sources in the northern part of the region which 
are marked in white in Figure 31 The northern- 
most of these is the source marked in yellow in 
Figure [H By examining the records associated 
with the affected region on the map, we can see 
that where a source saturates several neighboring 
pixels in the peak-up image, there is a bleed-over 
effect that reduces the response of the rest of the 
IRS chip (BCD image) in those rows, and may 
increase the response in neighboring rows. The 
result is broad, uneven dark or bright stripes in 
a 2D map extracted from the cube. The loca- 
tion varies with wavelength because as the PUI 
is scanned across the bright source in space, the 
artifact changes position on the BCD frame, i.e. 
changes wavelength. Simply masking all affected 
pixels at the record level results in holes in the 
maps at many wavelengths. Instead, we chose to 
mend the response. We selected all rows that in- 
tersected a group of saturated pixels on the peak- 
up image, with a margin of two rows on either 
side. These are the affected rows. We found the 
continuum level in adjacent rows, excluding major 
emission lines, and linearly interpolated across the 
affected rows. In the case of affected rows that cor- 
responded to continuum emission, we replaced the 
values with the interpolated values. In the case of 
affected rows that corresponded to atomic or PAH 
emission, we attempted to restore the base contin- 
uum level below the emission line by adding the in- 
terpolated value, and subtracting out the median 
value of the affected rows in the region between 
the SLl and SL2 chips, which is the bias level of 
the affected rows. The artifact is reduced, though 
not eliminated. The result is a complete cube, 
but with localized regions that have larger uncer- 
tainty in the SLl and SL2 modules (see the cyan 
and magenta regions marked in Figure [4]). For the 
purposes of presentation, the areas where the arti- 
fact 's effect is still distinct have been cropped out 
of each map. 

At some wavelengths in the SL cubes, a large 
region of somewhat elevated emission is still visi- 
ble after our corrections. This artifact appears to 
be another result of the saturation in the peak- 
up images. After a saturation event, the response 
of the chip is elevated for the remainder of the 
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AOR. The resulting artifact, which can be seen 
in maps at a variety of wavelengths (for exam- 
ple, 5.2, 6.3 and 10.7/^m), is localized in space. 
This low-level artifact is also generally reduced in 
continuum-subtracted images, because the effect 
is broad on the BCD frame, and thus broad in 
wavelength as well. 

2.3. Spectral Features and Line Fitting 

Hot dust emission is responsible for the over- 
all shape of the continuum in 30 Doradus and its 
strong rise to ward the red. U sing the fitting pack- 
age PAHFIT (jSmith Il2007ah . we fit the continuum 
with an assortment of four thermal functions from 
40 to 300K and a stellar blackbody of 5000K (the 
fit is insensitive to exact temperature as long as it 
peaks blueward of our shortest wavelength 5/im). 
The integrated spectrum of the entire region is 
shown below in Figure and astrophysically in- 
teresting spectra are described below in H3.3\ 

The broad PAH bands and unresolved atomic 
lines are blended in some parts of the spectrum. 
In order to decompose them, we fit the whole 
suite of emission lines at once. Detailed analy- 
sis of the dust features will be addressed in a fu- 
ture work. The feature of greatest potential con- 
cern is blending of [Ne II] A12.8/im with the 12.7/im 
PAH feature. The PAH contribution is much lower 
than the ato mic line over much of 30 Doradus. 
In particular, ISturm et al. I (|2000l) found that at 



the higher spectral resolution of ISO-SWS, the 
12.7^m PAH feature was undetectable relative to 
[Ne II]A12.8^m within the two /S'O-S' 1^5" apertures 
(located on the ridge) - the lowest ratio of PAH to 
[Ne II] of any H II region they studied. We used 
comparisons with high-resolution spectra (see Ap- 
pendix |B]) and tests in parts of 30 Doradus that 
should be completely PAH-free, like the core of 
R136, to determine that the joint fitting of those 
two features is robust in 30 Doradus and uncer- 
tainties are properly accounted for. 

2.4. MIPS SED cube 

We also made contemporaneous observations!! 
of 30 Doradus using the SED mo de of the MIPS 



(jRieke et all 120041: ILu et al.ll2008h . The raw data 
are responsivity-corrected, dark-subtracted, and 
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Fig. 4. — Three-color log-scale image generated 
from the IRS cube (see §2.ip with regions of in- 
terest overlaid. Black contour: a single radio 
contour to guide the eye (see Figure The 
emission line maps are continuum-subtracted, and 
each is plotted at a scale that brings out detail. 
Red: [S III]A33.4^m, in the LLl module. Green: 
[S IV]A10.5/xm, in the SLl module. Blue: PAH at 
6.2/im, in the SL2 module. The areas outlined in 
bold black are vulnerable to falling response (sat- 
uration) and strong fringing in the LLl module 
(see Appendix [B|) . Continuum fits in these areas 
should be regarded with caution, as should the fit- 
ted line strengths of long-wave emission lines. The 
sources marked in white saturated the peak-up 
camera and caused dark stripes and other artifacts 
in the SL modules elsewhere in the map. The cyan 
regions represent the total affected area for the 
SLl module, while the magenta regions represent 
the total affected area for the SL2 module. How- 
ever, the main effect of the saturation is a stripe 
across the map that changes with wavelength. 
There are generally no severe effects in the rest 
of the map for a given wavelength. The approxi- 
mate location of this stripe is marked by a series 
of bold lines, for several wavelengths correspond- 
ing to major emission lines in our spectra. Blue: 
[Ar n]A7.0/xm. Cyan: PAH at 6.2/zm. Magenta: 
Humphreys-7 A5.90/im. Green: [Ar III]A9.0/xm. 
Yellow: [S IV]A10.5Mm. Red: [Ne II]A12.8Mm. See 
^2.21 for more detailed discussion. 



illuminat ion-corrected using; t he MIPS germanium 
pipeline (jGordon et al.ll2007l ). with the difference 
that the illumination correction is derived from 
Galactic cirrus rather than zodiacal light. Wave- 
length correction is derived from observations of 
bright planetary nebula. Since 30 Doradus is so 
bright, we use the chopped "background" area 
simply as additional on-source observation, and 
assume that the thermal background is much less 
intense than the emission from the nebula itself. 
As with the IRS data, mapping the area with off- 
sets equal to half of the slit width, perpendicular 
to the slit direction, yields a fully sampled spec- 
tral cube, with resolution X/6X ~25. Even with 
such poor resolution, we detect [N III]A57.3/im 
with high significance over most of the map, and 
tentatively [O l]X63.1iim in the very center f ijS.ip . 
The shape of the dust continuum will be dis- 
cussed in a future publication on the dust content 
of 30 Doradus, but we note that the SED peaks 
shortward of 75/im, consistent with very high av- 
erage dust temperatures (Figure [S]). 

3. Results 

Figure [5] shows the integrated spectrum of 
30 Doradus, with dust and atomic features la- 
beled. The overall continuum shape is typical 
of hot small dust grains in H II regions. The 
MIPS/SED data clearly show that the spectral en- 
ergy distribution of the entire nebula peaks short- 
ward of 70^m, indicating a quite warm average 
dust temperature. A full discussion of dust prop- 
erties will be the subject of a subsequent publica- 
tion. We note that the aromatic features are well- 
detected but not extremely strong - comparable 
to Galactic compact H II regions or compact blue 
galaxies, and weaker than the integrated spectra 
of entire spiral galaxies or the diffuse emission of 
the Milky Way iGalhano et al.ll2008l : [Peeters et af 



show spectra of notable compact regions. 
3.1. Spatial distribution of emission 



[2002). The equivalent width of the 6.2/xm feature, 
averaged over this large aperture, is 0.61/im. In 
particular, the weak 17/im feat ure may indic ate 
a dearth of large neutral PAHs (jSmith Il2007al ) in 
this intensely ionized environment. The strongest 
spectral features are the ionic emission lines which 
are sensitive to physical conditions in the ionized 
gas. 

In the following sections we describe the spatial 
distribution of emission lines and their ratios, and 



Figures [6l ll4l show maps of ionic line emission 
in 30 Doradus. To first order, line emission fol- 
lows the diffuse emission pattern seen at other 
wavelengths, of a broad arc or "ridge" surround- 
ing an evacuated hole. The overall excitation of 
the region is immediately clear from the distri- 
bution of [Ar II] and [Ar III] - there is no de- 
tected [Ar II] in the evacuated hole. With an 
ionization potential of 15.7eV, the photoioniza- 
tion cross section and charge transfer with hydro- 
gen e nsures very little Ar^ within an H II region 
(e.g. Sofia fc Jenkins! Il998l ). (Ionization poten- 
tials are listed in Tabled] for reference.) We note 
that the Ar"'"^ recombination rate is highly un- 
certain, so absolute calibration of [Ar III]/[Ar II] 
has large systematic uncertainty ( Morisset et al.l 
2004 IStasihska k Schaereilll997n . The similarity 
of overall morphology amongst the Ar, Ne, and 
S lines suggests that any variation of excitation 
across the nebula will be correlated with the struc- 
ture seen in continuum from radio to infrared, al- 
though as we will explore thoroughly below, the 
shape of the hole varies somewhat, and there are 
some distinct locations of high excitation. The 
[Si II] distribution appears somewhat displaced 
away from R136 and the cavity, as would be ex- 
pected since the line is strong in PDRs as well as 
ionized gas. 

Lower ionization species have tentative detec- 
tions in a few places, but never very strongly 
or widely distributed. Of particular interest is 
a peak near 5''38"44.9" -69°Q5'13" (source "0", 
i )3.3p . at which we also detect [N III]A57.3/im 
and [O I]A63.1/im in MIPS SED observations (see 
the northeastern region marked in Figure llip . 
The spectra in this part of the map are affected 
by fringing in the LLl module (>20/im), so a 
[Fe II] A26.0/im detection in the same region is con- 
sidered tentative. The ratios of [Ne III]/[Ne II] and 
[S IV] /[S III] (and the MIR continuum) are high 
in this general area, but this region in the center 
of the ridge is complex, and the higher ionization 
species do not peak in exactly the same location 
as these lower ionization species. In fact, the peak 
of low-ionization species falls between two peaks 
of the [Ne III]/[Ne II] map. 

Molecular hydrogen emission is detected in the 
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Table 1 
Ionization Potentials. 
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II 


III 


IV 


o 


13.62 


35.12 


54.94 


77.41 


Nc 


21.56 


40.96 


63.45 


97.11 
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10.36 


23.33 


34.83 


47.30 


Ar 


15.76 


27.63 


40.74 


59.81 



Note. — Units arc cV. 
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Fig. 5. — The spectrum of the entire 30 Doradus 
nebula (a square about 4' or 60 pc on a side cen- 
tered on 05'*38"44^ -69°05'28"). The observed 
spectra from IRS and MIPS/SED, and infrared 
photometric points, are plotted with symbols, 
while the fit we made to the IRS data is displayed 
with a solid line. 



central molecular ridge and marginally detected in 
several other locations. Only the S(3) 9.67/im line 
is reliably detected over much of the map. Tenta- 
tive detections in our low-resolution spectra of the 
S(0) 28/im line can be ruled out by examining the 
handful of high-resolution spectra in 30 Doradus. 
Detection of the S(2) 12. 28/im line independently 
of the neighboring Humphreys-a A12.37/im line is 
not very reliable with low-resolution spectra as de- 
scribed in more detail below, and from the high- 
resolution spectra we find that the molecular line 
is weaker than the atomic line everywhere, and 
less than 15% over most of the central part of the 
map. 

Higher excitation species [O IV] and [Ne V], 
which have strong lines in the IRS bandpass, are 
not detected in our maps. As discussed in tJ5.21 
this favors photoionization over shocks as the dom- 
inant physical process in the region. 
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G4.a00 54.750 G4700 S4.650 94 600 34 550 



Fig. 6. — through Figure 1111 Fitted atomic hue maps. AU plots are on the same log scale, ranging from 
10^^ to 10^''' Wm~^sr~^ . Note that four of the fainter maps have been multiplied by 10, and the brighter 
[O III] map has been divided by 100. The length of the black bar is 1" or 15 pc. White contours: a single 
level of 3cm radio emission, to guide the eye (see Figure [2|). 
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Fig. 7. — Fitted line maps (see caption. Figure [H]). 
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Fig. 8. — Fitted atomic line maps (see caption. Figure [6]). 
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Fig. 9. — Fitted atomic line maps (see caption, Figure [H]). 



[S III] 18.7 urn 




[Ne III] 15.5^m 
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Fig. 10. — Fitted atomic line maps (see caption, Figure[6|). The [S III]A18.7/im map is marked with regions 
of interest listed in Table [2l 
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Fig. 11. — Fitted atomic line maps (see caption. Figure [H]). The black circles on the [O III]A88.3/im map 
mark the detected peaks of [O I]A63.1/im and [N III]A57.3/im. 
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PAH 6.2 nm 
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Fig. 12. — through Figure [HI Fitted PAH feature maps. For details, see caption, Figure [SI 
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Fig. 13. — Fitted PAH feature maps (see caption. Figure [5]). 
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Fig. 14. — Fitted PAH feature maps (see caption. Figure [H]). 
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Fig. 15. — Line ratios on a linear scale. White contours: a single level of Scm radio emission, to guide the eye 
(see Figure[2]). The white region in the middle of the argon ratio plot indicates nondetection of [Ar II]A7.0/.tm 
and a very high value of the ratio. Small masked regions in the lower left of the plots are affected by artifacts 
(see Fig. g]). 
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3.2. Line Ratios 

Figure [15] shows the three Hne ratio maps 
[Ar III]A9.0/im/ [Ar II]A7.0/xm, [Ne III]A15.5^m/ 
[Ne II]A12.8Aim, and [S IV]A10.5/zm/ [S III]A33.Vni. 
These species all have ionization potentials above 
13.6eV, and thus are sensitive to the shape or 
hardness of the ionizing spectrum (we will call 
this Trad J because it is often parameterized by the 
temperature of the best-fitting black-body over 
A <912A) and the intensity of ionizing radiation. 
The latter is usually quantified as the dimension- 
less ionization parameter 



U : 



1 



^ Jl3.6eV 



chv 



or ratio of ionizing photon density to atom density. 
Any single line ratio cannot distinguish between 
elevated U or Trad, but because the spacing of the 
ionization potentials differs with atom, different 
ratios have different U and Trad dependencies, and 
measuring two ratios can break the degeneracy. 
(This has been discussed b y many authors, see 



especially the discussioii s in Moris set et al 



2004 : 



20061 



iMartfn- Hernandez et al. [2002; .Dopita et al 
,and gSDand Fig HD] below.). 

All ratios increase with [/, but locations where 
the ratios are not well correlated can indicate 
changes in the hardness of the ionizing field. In 
30 Doradus the ratios are very well correlated (cor- 
relation coefficient of 0.7 between the neon ratio 
and argon ratio map, 0.8 between the neon and 
sulphur ratios, and 0.7 between the argon and sul- 
phur ratios). The Ne and S ratios are telling us 
similar things over much of the region. For ex- 
ample, there is high excitation in a "hot spot" 
between the two lobes of the ridge (source "A", 
i)3.3p . and in the region to the south of the bub- 
ble. There is a low excitation ridge to the north 
(source "F", §3.3p . The location noted above for 
low argon emission (source "G", ^3.3p is signifi- 
cantly low in all three line ratios. 

Interestingly, the Ne and S line ratios show 
quite different behavior on the eastern edge of 
the bubble, where [Ne III]/[Ne II] is high but 
[S IV] /[S III] is low (the red region in Figure [H] 
below, centered on source "E", tj3.3p . The S ratio 
is more affected by extinction since [S IV]A10.5/im 
is in the silicate dust absorption feature, so that 
could be a region of very high extinction or a re- 
gion of particularly hard ionizing radiation. In 



this case, the sulphur ratio does not vary dra- 
matically throughout the bubble, and it is the 
neon ratio which is higher at the western end, 
which argues for harder ionizing radiation. An 
extinction effect would require simultaneous in- 
crease in the ionization parameter or strength of 
the ionizing field and increased extinction. The 
specific extinction difference (A([S IV]A10.5/im)- 
A([S III]A18.7//m))/A(2.1/^m) equals 0.42±0.03 
(uncertainty reflects differences between published 
extinct i on cu rves; we used that in PAHFIT, 
ISmith I l2007a[ ). whereas the effect of extinction 
on the neon ratio is small (A([Ne III]A15.5/im)- 
A([Ne II]A12.8/im))/A(2.1/xm) = 0.0±0.07. Re- 
producing the observed line ratios without chang- 
ing the hardness of the ionizing field would require 
A(2/im)~0.75, or Ay ~6. Such significant extinc- 
tion seems unlikely given the overall relatively low 
extinction in the region, and is not detected in our 
extinction maps (t j4.ip . 

One question that is possible to investigate with 
this dataset is whether the optically known hot 
stars in 30 Doradus can completely account for the 
ionization structure, or whether embedded star 
formation affects gas excitation measured in the 
infrared that otherwise might be concealed by at 
other wavelengths. On Figure [16] we plot show 
the location of Wolf-Rayet stars and early O stars 
(there are discrepancies between spectral types de- 
termined by different authors, and our intent is 
to show the most massive stars, not provide the 
most precise spectral typing possible). There is 
a cluster of WR stars between the two lobes of 
the ridge, coincident with very high excitation gas. 
There is not a dramatic increase of excitation cen- 
tered on R136, although the excitation is generally 
high in the area. Figure [T6] also show s the 3 .5/im- 
excess sources of iMaercker fc Burtonl (|2005l ). rep- 
resenting a crude selection of possible protostel- 
lar candidates. Other kinds of sources, including 
some of the previously identified WR stars, can 
also display infrared excess. Nevertheless, this se- 
lection of sources does trace what is known from 
more pre cise studies w ith incomplete spatial cover- 
age (e.g. lRubio et al.l ll998: Brandncr et al. 200ll), 
that the young embedded sources in this region 
are concentrated along the IR-bright ridge. The 
protostellar candidates shown here do not show 
any particularly striking correlation with regions 
of high excitation - in fact the southern part of 
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the region has high excitation and few protostel- 
lar candidates. We conclude that the ionization 
structure in 30 Doradus is primarily determined 
by the (optically) known hot stars. 

3.3. Notable Regions 

Figure [T7] shows the spectra of small regions 
which stand out in the continuum (Fig[T]), in the 
feature maps (Figures [6l through [T4)) . or in the line 
ratio maps (Fig [T5)) . They have been chosen to 
illustrate the range of excitation, extinction, and 
continuum shape in the 30 Doradus nebula. These 
regions of interest were marked on Figure [21 and 
are listed for reference in Tabled] 

Source "0" is the brightest point on the ridge 
in most tracers, including infrared continuum, 
many lines, and centimeter continuum. We de- 
tect [N III]A57.3Aim and [O l]X65.1^J.m in MIPS 
SED observations at this location, and tentatively 
[Fe II]A26.0/im, although this small region of the 
map is affected by residual fringing and saturation 
in the LLl module, raising the uncertainty of the 
[Fe II]A26.0//m detection. 

Source "A" is a diffuse region between the two 
lobes of the molecular ridge. We describe it as a 
"hot spot" because it has locally elevated excita- 
tion seen in the line ratios (FigllSp. and when ana- 
lyzed with photoionization models, is best fit with 
a hotter or harder radiation field (Fig [^11 > j5.ip . 
There are known hot stars in the vicinity, but it 
is not completely clear if they are responsible for 
the locally higher excitation in the H II region. 

Source "B" is a prominent point source from 
optical to MIR wavelengths, and has bee n spec- 
trosco pically classified as an M supergiant I Parker 
( I993L ; see also refs therein). Their quoted ab- 
solute V magnitude appears to place the star at 
twice the distance to the LMC, but other authors 
note that it is likely in the foreground. It seems 
unlikely in any case that it is located in the bubble 
or affecting the ionization structure. Because this 
source exhibits silicate emission at 9.8 and 18 /im 
rather than absorption, PAHFIT cannot be ap- 
plied to it using the parameters used on the rest 
of the spectral cube. It is the only source in Fig- 
ure [17] without the fitted spectrum overlaid on the 
data points and errorbars. 

Source "C" is the R136 cluster core. 

Source "D" is a prominent mid-infrared point 



source outside the main cluster and bubble. The 
surrounding ionized gas has quite high atomic line 
ratios and thus appears quite highly excited, but 
we were unable to find a hot star cataloged in the 
literature which might be responsible for this local 
excitation. This region also shows relatively high 
extinction compared to the rest of 30 Doradus. 
Extinction can raise the sulphur line ratio but has 
relatively small effect on the neon line ratio, so in 
our opinion the extinction is not responsible for 
the apparently high excitation. 

Source "E" is a known WN6 star R145, fairly 
isolated from other catalogue hot stars at the east- 
ern end of the bubble. The high fitted radiation 
temperatures in the neighborhood of this source 
( i|5.ip are consistent with excitation by the harder 
expected radiation field from a Wolf-Rayet star. 

Source "F" is an east- west extended region that 
appears to be a trough or low point in the ex- 
citation and line ratios. The electron density is 
not remarkable, neither does there appear to be 
any molecular gas traced by CO emission in the 
vicinity. It is possible that the feature is in the 
foreground of this three-dimensional nebula. 

Source "G" is a region associated with two 
pea ks in the CO emission [JGB98] 30 Dor-12 and 
13 ( Johansson et al. 19981 ). just to the west of the 
main ridge. There is a peak in both argon line 
maps at this location, and the region is a low 
point in all three excitation ratios [Ar III]/[Ar II], 
[Ne III]/[Ne II] and [S IV] /[S III]. There are in fact 
several discrete parsec-sized regions of low excita- 
tion in the vicinity There is little 3cm continuum 
at those locations, which lie almost exactly on the 
opposite side of a molecular cloud from R136. One 
spot coincides with a knot in the infrared diffus e 
continuum (#29) noted bv iHvland et all (|l992l ). 



There is a 3.5jnm-excess protostellar candidate 
( Maercker fc BurtonI l2005l and Fig. [H]) nearby, 
but not coincident with the low-excitation spots. 
High resolution and sensitivity molecular observa- 
tions may reveal dense, starless molecular clumps 
in this part of 30 Doradus, (self )-shielded from the 
intense radiation. 
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Table 2 
Notable Sources 



label 


RA 


Doc 


description 


"0" 


5h38m45s 


-69d05ml3s 


brightest point on ridge 


"A" 


5h38m42s 


-69d05m05s 


hot spots between molecular clouds 


"B" 


5h38m49s 


-69d05m32s 


Parker (1993) #1445 MI star 


"C" 


5h38m42s 


-69d06m00s 


R136 


"D" 


5h38m56s 


-69d04ml8s 


example of high extinction 


"E" 


5h38m57s 


-69d06m07s 


isolated WN6 star R145 


"F" 


5h38m57s 


-69d03m42s 


low-excitation trough 


"G" 


5h38m28s 


-69d06m30s 


SW low-excitation region 



Distribution of Matter: Extinction and 




Fig. 16. — IRAC 8^m image (log scale) with two 
types of relevant point sources: Circles show the 
most energetic stars in the reg ion - Wolf-Rayet 



stars a nd early O s tars, from iBrevsacher et al 
(|l999t ): IParkeil (|l993l ). and X's show protostellar 
ca ndidates determined from S/^m excess emission 
bv lMaercker Sz BurtonI ( 20051 ). White contours: a 
single level of 3cm radio emission, to guide the eye 
(see Figure 



and 
of 



4.1. Extinction 

Previous estimates of t he extinction in 30 Doradus 
found Av=l.l±0.1mag (|Dickel et al.l lit 
refs therein), with two possible locations 
higher extinction at 05^ 38'"39" -69°0r3 " and 
05''38™32" -69°06'22". iRosa fc Mathi^ (|l987l) 
found extinction corrections at H/3 of 0.5-1.0 
(Av'=0.4-0.8) at ten locations in the outer parts 
of the region we are studying. We make three 
estimates of extinction in the region, two directly 
from our dataset. All methods suffer from system- 
atic uncertainty and modest signal-to-noise, but 
we can be confident of higher extinction in the 
regions that all three methods agree. 

First we use the ratio of centimeter contin- 
uum to Ha emission to derive a map of extinction 
in 30 Doradus, follow ing the procedure used in 



Lazendic et al.l (|2003l ) and using their centimeter 
data that they kindly provided (magenta contours, 
in Fig. [T5)) . We assume that all of the centimeter 
continuum emission in the 30 Doradus region is 
thermal. Single dish measurements estimated that 
the nonther mal component con tributes less than 
2% at 6cm ( Shaver e t al. 19831 and refs therein). 
Lazendic et al . (2003) identified two possible su- 
pernova remnants from comparison of cm synthe- 
sis images and optical recombination lines, but 
their own analysis and subsequent foUowup with 
optical and Xray imaging and optical spectroscopy 
suggests that these arc merely extinguished H II 
regions (Chu et al. 2004). The relationship be- 
tween thermal bremsstrahlung centimeter contin- 
uum and hydrogen recombination line emission 
depends weakly on th e ele ctron temp e ratur e, but 
Rosa fc Mat"hi3 (|l987l ) and lPeck et all (|l997t ) both 
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Fig. 17. — Spectra of the seven regions described 
in H3.31 with errorbars. For comparison, the spec- 
trum of the entire region (Fig [5]) is repeated at the 
top, with feature labels. 



found that Te variations are small in 30 Doradus 
(< ±300K) , so it is unlikely that the calculated ex- 
tinction variations are actually misinterpreted Te 
variations. 

We derived two more maps of extinction di- 
rectly from the IRS spectral cube. Amorphous 
silicate dust is responsible for two broad bands of 
absorption at 9.7 and 18^m. The PAHFIT pack- 
age reports the fitted optical depth at 9.7 /im. The 
resulting map for 30 Doradus is shown in Figure 
[TSl in green contours. The fit to the absorption 
feature is sensitive to noise in the spectrum. We 
cropped away southern parts of the map that were 
clearly artifacts, and any part of the map where 
the signal-to-noise was less than 2. Finally, we 
median-smoothed the resulting map with a win- 
dow of 3 pixels. The result is a sparse map show- 
ing the regions of comparatively reliable elevated 
silicate absorption. 

The ratio of Hydrogen recombination line 
strengths can also be used to calculate extinc- 
tion, most commonly by assuming an intrinsic 
(unextincted) ratio from Case B recombination. 
Our dataset includes Humphreys-a (7-6) and 
Humphreys-7 (9-6). Due to modest signal-to- 
noise and to a systematic tendency of PAHFIT to 
slightly underestimate the continuum level near 
these wavelengths, the absolute value of extinction 
calculated from this ratio has a large systematic 
uncertainty. It is also difficult to resolve the weak 
Humphreys-a line from H2 S(2) 12.28/im in low- 
resolution spectra. We quantified the amount of 
potential contamination to the Humphreys-a line 
strength in several ways: We fit our low-resolution 
spectra with and without the molecular hydrogen 
line (the central wavelengths of both lines are 
very tightly constrained by PAHFIT), and found 
that the Humphreys-a line strength was only de- 
creased by 15% (that flux was attributed to H2) 
when both lines were included. We also fit the 
high-resolution GTO spectra, in which the lines 
are easily separated and reliably measured, and 
found that over most of our area, H2 was less 
than 20% of Humphreys-a. We examined fits to 
both the high and low-resolution spectra in the 
few high-resolution apertures where the H2 and 
Humphreys-a strengths are comparable (H2 never 
exceeds Humphreys-a in strength) and found that 
the low-resolution Humphreys-a strength used in 
our analysis was at most overestimated by 40% 
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by the presence of the H2 Hne on the wing of the 
Humphreys-a Hne. While some level of contami- 
nation of Hmnphreys-Qf by H2 S(2) may be present 
in our maps, we do not expect that to change any 
of our conclusions including the regions of high 
extinction identified in Figure [TSl 

Further systematic uncertainty arises in the 
choice of extinction curve, since the lines lie 
on the wings of the silicate absorption features, 
which vary in amplitude and width in different 
studies. However, all extinction curves that we 
considered show greater relative extinction at 
Humphreys-a A12.37/im than at Humphreys-7 
AS .QO^m. We e xamined the extinction curves 



of IChiar fc Tielens (,2006) for the Galactic Cen- 
ter and the ISM, and the average LMC extinc- 
tion curv^ base d on the carbonaceou s- silicat e 
grain model of IWeingartner fc Draind (|200l[ ). 
All of the extinction curves roughly agree in 
the vicinity of the Humphreys-a emission line. 
Near the Humphrcys-7 line, the Chiar fc Tielens 
curves agree well with t he near-infrared ex t inctio n 
for the ISM found by llndebetouw et all (|2005h . 
Thus, extinction correlates with the flux ratio of 
Humphreys-7 A5.90/J,m to Humphreys-a A12.37/im, 
and we can still find the relative level of extinction 
across 30 Doradus from the Humphreys-a/Humphreys-7 
ratio, even though the absolute normalization is 
uncertain. 

We calculated the signal-to-noise of the ratio 
based on the RMS variance of the individual maps 
of emission line strength, and removed all re- 
gions of the map where the signal-to-noise was 
less than 2. We cropped away the parts of the 
map where an additive correction was made to 
the SL spectra (see iiA.0.3|) . because this subset 
of the map encloses some spectra which have an 
unphysical plunge at the red edge of the SL2 mod- 
ule. The fit to those spectra generally underesti- 
mates the continuum near the Humphreys-7 line, 
and thus overestimates the ratio of Humphreys-7 
to Humphreys-a. We also cropped away the parts 
of the map where the reduced oi the fit was 
greater than 5. Finally, even after applying these 
filters, it was necessary to trim away some regions 
from the low-signal edge of the map where the 
noise levels were still untrustworthy. The ratio 
map shown in grayscale in Figure [18] is now mainly 



limited to those areas of 30 Doradus where we ob- 
tained both good detections of these two faint lines 
and a good fit. 

Despite the level of noise in the maps of ex- 
tinction, particularly the map derived from the 
Humphreys-a and Humphreys-7 lines, there are 
three areas where all three methods agree on es- 
pecially high extinction. These areas have been 
marked on Figurc fTSl Generally, there is not much 
extinction by dust in 30 Doradus. As discussed 
above, extinction has a small effect on the Neon 
excitation ratio, but raises the [S IV] /[S III] (due 
to the shape of the extinction curve) , which in turn 
will get interpreted as higher Trad in our photoion- 
ization models. The regions of interestingly high 
Trad do not turn out to correspond to regions of 
clearly high extinction, but the effect should be 
kept in mind in interpreting the data. 

4.2. Electron Density 

Figure [19] shows the electron density calcu- 
lated from [S III]A18.7^m/[S III]A33.4/im. The 
[S III]A18.7/im map was convolved to the lower 
resolution of the [S III]A33.4/im map, and the 
line ratio conv erted to elect r on de nsity using the 
conversion in iDudik et all (|2007[ ) at Te=10^K 



(their Figure 9 and section 6: Those authors 
calculated the line ratio as a function of den- 
sity and temperature for a five level atom us- 
ing th e collision strengths from iTaval fc Gupta 
(I1999) and ra diative transition probabilities from 
Mendoza fc ZeippenI (jl98^ ). The calculated ra- 
tio is not a strong function of temperature. The 
ridge is prominent in the [S III] ratio or Ug map, 
and in fact the Ug map is quite similar to the 3cm 
continuum in morphology. The density is elevated 
to the south of the bubble and R136, in a region 
of relatively high excitation. 

4.3. Abundances 

Variations in elemental abundances can in 
principle change the diagnostic ratios that we 
are using to measure physical conditions in 
the 30 Doradus nebula. Previous studies have 
found very small internal abundance variations 
in 30 Do r adus and other giant H II regions. 



^Available at: www.astro.princeton.edu/ draine/dust/dustmix.html 



Peck et al.l (|l997f l found no significant variations 
on 15" scales in the He abundance Y'^ = He'^ /H~^ 
— 0.13±0.02 measured from radio recombination 
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Fig. 18. — Map of extinction in 30 Do- 
radus. Grayscale linear image: the IRS 
Humphreys-7/Huniplireys-a ratio, which is pro- 
portional to extinction. Dark areas of the 
map indicate high extinction. The Humphreys 
lines are weak, and this map has been severely 
cropped to remove areas of the worst signal-to- 
noise. Green contours: the optical depth of sili- 
cate absorption from the IRS spectral map, r — 
0.10,0.27,0.43,0.60. Areas of low signal-to-noise 
have been masked out of this map as well. Ma- 
genta contours: The ratio of the 3cm continuum 



(jLazendic et all l2003h to Ha emission (MCELS, 



http://www.ctio.noao.edu/~mcels/), also propor- 
tional to attenuation. The northernmost ma- 
genta source corresponds with the molecular cloud 
marked in yellow in Figure [TJ White contours: 
a single level of 3cm radio emission, to guide the 
eye (see Figure [2]). The areas where all three maps 
tend to agree on high extinction have been marked 
in black. The northeasternmost corresponds with 
source "D" in Figure [2] and §3.31 



lines. Rosa fc Mathid ( 1987 ) found little variation 
between several positions measured with optical 
spectr oscopy. Of m ore direct relevance to this 



work, iLebouteiller et al.l ( 20071 ) found less than 
0.01 dex dispersion of Ne/H, S/H, and Ar/H, 
using the high resolution GTO Spitzer spectra 
mentioned above. 

Small scale abundance variations in H II re- 



gions including 30 Doradus (jTsamis fc Peauignot 
2OO5I ) have been proposed to explain discrepancies 
between optical and infrared abundance determi- 
nations. We might hope to detect sub-parsec scale 
abundance variations in this dataset by using three 
infrared line ratios to solve simultaneously for f7. 
Trad, and Z . In practice, we find that the excita- 
tion variations can be adequately explained with- 
out abundance variations, which would be seen as 
residuals in our fitting of V and T^ad- The signal- 
to-noise especially in our argon line ratio is insuffi- 
cient to detect ab undance variations at the 0.1 dex 
level predicted bv lTsamis^fc Peauignot (2005} ). 

Overall, abund ances in 30 Do r adus are not 
particularly low: Rosa fc Mathid ( 1987 ) found 
Ne/0, S/0, and Ar/0 ratios close to solar. 
Lebouteiller et all (|2007l ) found Nc/H and Ar/H 
of 12-Hlog(X/H)=7.76±0.02 and 6.32±0.06, re- 
spectiv ely, within 0.25 dex of the range of sola r 
values (|Lodderdl200"7l l2003t lAsplund et al.ll2005^ . 
They found S/H of 6.77±0.03, only about 0.3 dex 
sub-solar. We also find that on average, half-solar 
abundances result in modestly better agreement 
between the three infrared line ratios in our pho- 
toionization models than solar or 0.1-solar models. 
In practice, systematic effects such as the argon 
recombination rate have a larger effect on this 



agreement than the abundances (jMorisset et al 
2004HStasinska fc Schae"mlll997l ). 



5. Distribution of Radiation and Gas: Ex- 
citation 

5.1. Photoionization models 

As mentioned above, if photoionization is as- 
sumed to be the dominant physical process, the 
ionic line ratios depend on U, Trad, and to a lesser 
degree metallicity. If one assumes constant abun- 
dances and hardness or Trad, then the line ratio 
maps in Figure [15] are maps of the ionization pa- 
rameter U, varying from log U ~-3 to log U ~-1.5. 

A somewhat more sophisticated analysis is to 
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solve for Trad and U simultaneously. We prepared 
a grid of photoioni zation models using Cloudy 



(jFerland et al.l 119981 ) as a 0-dimensional tool to 



solve for the ionization structure and line emissiv- 
ities given a specified radiation field and ionization 
parameter (we used the output of the first zone 
in each simulation). We exp lored different in- 
put s pectra, in cluding ATL AS ("Castelh fc Kuruc3 
20041 ). Tlusty (|Hubenv fc L anz 1995), CoStar OB 
and Wolf-Rayet atmosp heres (jSchaerer fc de Koter 
19971: ISmith et al.ll2002f) . and black-bodies. In the 
end we used a grid of Tlusty atmospheres cal- 
culated at half-solar metallicity, extrapolated to 
hotter effective temperatures using black-body at- 
mospheres to set the functional dependence of 
line ratios on effective temperature, and normal- 
izing the ratios to those of the hottest Tlusty 
models. (We found that adopting different stel- 
lar atmospheres primarily changes the line ratios 
by a constant multiplicative factor, and has very 
little effect on the functional dependence of the 
line ratios on effective temperature and ioniza- 
tion parameter.) We were not able to find freely 
available grids of more modern atmospheres (e.g. 
WMBASIC, CoStar) at subsolar metallicity, but 
we performed careful comparisons of our modeling 
at solar metallicity to understand the systematic 
effects. If nebular abundances are set at solar 
levels in the photoionization models, and stellar 
atmospheres calculated at solar abundance are 
used for self-consistency, the derived ionization 
parameter U decreases systematically by about 
0.1 dex, and the derived radiation temperature 
increases systematically by 10%. At solar metal- 
licity, WMBASIC atmospheres result in about a 
factor of two lower [Ne III]/[Ne II] ratio and a 
50% lower [S IV]/[S III]ratio. That would increase 
the derived ionization parameter in 30 Doradus 
systematically by ^0.15 dex, and lower the de- 
rived radiation temperature by 5-10%. These 
effects have been explained in d e tail by other au- 
thors, especially iMorisset et al.l (2004). None of 
the changes in stellar atmospheres or metallic- 
ity that we explored would result in qualitative 
changes in our conclusions, merely small system- 
atic shifts in derived parameters. We also varied 
the dust prescription in terms of abundance and 
grain size. Neither had a strong effect on the line 
ratios as a function of T^ad and provided that 
U was calculated locally, i.e. from integrating the 



diffuse ionizing field in the simulation above 1 Ry, 
rather than assuming some geometry-dependent 
expression such as Ni^/c^-Kr^n,,- 

Figure [^D] shows typical behavior - each mea- 
sured ratio defines a curve in [/ — T^ad space, but 
the lines have different slopes because the ioniza- 
tion potentials are spaced differently for the differ- 
ent atoms. This is equiv alent to the nea rly parallel 
arro ws for U and Tmd inlMorisset et al.l (12004,1 fsee 



also Martin-Hernandez et al.l l2002t iDopita et al 
20061 ). Two measured ratios can better constrain 
U and Trad- We will focus on the Ne and S ratios 
since [Ar II] is not detected over large parts of the 
H II region. 

Figure [H] shows the fitted U and Trad across 
30 Doradus. The ridge is a region of high ioniza- 
tion parameter, as previously noted from simple 
examination of the line ratios. The region around 
the relatively isolated source at 5''38"*56.5* - 
69° 04' 17" to the north of the ridge is also highly 
excited, probably due to the local effects of that 
star (source "D" , Figure O §3.3p . We note two 
regions of high U in particular. The most promi- 
nent is between the two parts of the ridge (just 
south of source "A" , ij3.3p . where young stars may 
be locally ionizing the gas and beginning to lower 
the density and disperse the ridge (see Figure [TBI 
for the location of protostars and the most ener- 
getic optically identified stars). Alternately, the 
density is simply lower there and ionizing radi- 
ation can more easily leak out from the bubble 
region around R136. The northern side of this 
"hot spot" shows evidence for hardening of the 
radiation field, which could result of the radia- 
tion originates in the bubble or on its rim, and 
is propagating northward, and being absorbed by 
gas and dust (both of which will harden the field) . 
The high degree of porosity and mixing between 
molecular and ionize d material in 30 Dora dus is 
well known (see e.g. IPoglitschet aDll995h . On 
the south side of the bubble is another region of 
high ionization parameter - this region also shows 
high electron density in the [S III] ratio, so the 
field must be locally strengthened, perhaps by the 
WR stars known in that region (Fig [T6|) . Par- 
ticularly interesting is the region on the eastern 
side of the bubble, which shows up in the fitted 
parameter maps as high Trad, but not particularly 
high U. There is a single catalogued Wolf-Rayet 
star (R145, WN6 type, source "E", gSS]) in the 
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center of that area, which may be energizing the 
eastern end of the bubble. 

The radial variation of physical conditions in 
30 Doradus, as a function of distance from R136, is 
of particular interest, to determine whether feed- 
back in the nebula is dominated by that cluster 
core, or whether individual hot stars scattered 
throughout the region are equally important. We 
have already seen that the latter is true, for ex- 
ample the hot eastern end of the bubble appar- 
ently excited by the WN6 star source "E" . Within 
several tens of parsecs from R136, however, there 
does appear to be a global effect. The top panel 
of Figure [22] shows the radial dependence of the 
ionic line ratios and fitted T^ad and U . The bub- 
ble at ~10pc radius is clearly evident in lowered 
Trad and raised U . Outside the bubble there is 
a modest gradient in the hardness of the fitted 
radiation field, but not much gradient in the ion- 
ization parameter U . This is sensible in a region 
which has the hottest stars in the center, but some 
ionizing sources distributed more widely. As R136 
is apparently driving an off-center "blister" type 
H II region, it is interesting to separate the radial 
dependencies of the ionized gas parameters in the 
eastern and western direction. The second panel of 
Figure [22] shows this comparison, showing clearly 
the bubble wall at ^^Spc distance in the east, and 
the hot bubble between R136 and the '^25pc dis- 
tant bubble wall in the west. Unfortunately the 
effects of saturation make the measurements less 
reliable beyond the bubble wall in the west (see 



2007 ) and the stronger shocks in SNR (jNeufeld et al 



20071 ) show similar relatively low-excitation MIR 



5.2. Shock Models 

Protostellar shocks can also produce MIR fine- 
structure emission, but only if the excitation 
(shock velocity) is high en ough. [Ne II] requires 
J-sho cks with v>60km s~^ (jllollenbach fc McKed 
Il989l) and [Ne III] velocities in exce ss of 100km s~^ 
(iMolinari fc Noriega-Crespoll2002l ). iLefloch et cd 



poos') detect [Ne II]A12.8Atm and [Ne III]A15.5/im 
in HH 2, but only from the highest excitation 
working surface. Similar results are being found 
with Spitzer, detecting [Ne II]A12.8/im but not 
[Ne III]A15.5/;im, and [S III]A33.4Mm but not 
[SjV]A10.5/zm i n HH46/47 (iNoriega-Crespo et al ' 



[2004"), Cep E (iNoriega-Crespo et al 
HH7-11 and 54 (jNeufeld et al.l l2006[l 



2004 ). 



and 
Even the 

bow shock near a runaway 09.5 star ([France et al 



emission. 

The (few) MIR observations of regions con- 
taining both shocks and photoionization tend to 
show photoionization dominant. In [Ne II] and 
[S IV] maps of the massive star formation region 
W51 IR S2, high spectral resolution mapping with 
TEXES (|Lacv et al.[[2007[ ) allows the high velocity 
emission to be separated from that at the velocity 
of the molecular cloud. The authors do not note 
dramatically different line ratios in the high veloc- 
ity emission from that at the systemic velocity, but 
rather confirm earlier ground-based observations 
of MIR fine structure lines consi stent with pho- 



toion iz ation by lat e 0-ty pe stars (jOkamoto et al 
20011 ). [Lacv et al.l ( 20071 ) interpret their observa- 



tions as a neutral jet emerging from the molecular 
cloud and being subsequently photoionized. 

More recent work with S pitzer in the Galactic 
center ([Simpson et al.[ [20071 ) cannot produce the 
observed O"*"^ abundance with any reasonable pho- 
toionization model, even including hot supergiant 
atmospheres and 10^~^K blackbodies representing 
diffuse X-ray emission. The proposed alternative 
of ~100km shocks can match the observed 
[O IV] line emission. The authors do not require 
shock excitation to explain their Nc and S line ra- 
tios, but do note that the highest excitation gas 
is found between the two stellar clusters, which is 
indicative of shock excitation. 

We ran a se t of shock model s using John Ray- 



mond's code (Ravinond 



19791 : [Cox fc Ravmond 



1985[ : [Hartigan et al.l[l987l ) with a range of veloc- 
ities and pre-shock densities. The ratios of rele- 
vant fine-structure lines are shown in Figure 1231 
along with the observed ranges in 30 Doradus. To 
first order, all three ratios can be matched with 
shocks on the order of a few hundred km s~^ 
in quite low-density gas {n < Icm""^). How- 
ever, we have significant nondetections of emis- 
sion from the higher ionization species [O IV] 
and [Ne V], with 3-a upper limits less than 10^^ 
relative to [S III]A18.7/im across the map. The 
predicted ratios for those two species relative to 
[S III] are also calculated, and shocks of this order 
should produce easily detectable [O IV]. Shocks 
can a lso produce [Ne V] (e .g. in supernova rem- 
nants Rakowski et al 1 [2007t ). but it is more usually 
used as an indicator of very hard radiation from 
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AGN or WR stars (|Abel fc Satvapall l2008l and 
Table [1]). Thus we conclude that in 30 Doradus as 
elsewhere, although shocks certainly exist, the ion- 
ization structure of the gas is dominated by pho- 
toionization and not shock activity. Furthermore, 
while the effects of hard radiation from WR stars 
affects the ionization balance in some parts of the 
nebula, the radiation field on ~0.5pc scales is not 
apparently hard enough for very high ionization 
species to be very important. 

It is particularly interesting to consider the 
"hot" region on the eastern side of the bubble 
where the photoionization models are driven to 
high Teff and moderate U by low [S IV] /[S III] and 
moderate [Ne III]/[Ne II]. That low [S IV] /[S III] 
ratio makes this part of 30 Doradus more consis- 
tent with shock models than other parts, although 
the nondetections of [O IV] and [Ne V] are still 
problematic. We suggest that of anywhere in the 
nebula, the eastern end of the bubble may have 
the highest likelihood of excitation by shocks, as 
would be expected from winds of the WN6 star 
located there hitting the denser sides of an evacu- 
ated bubble. 
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Fig. 19. — Ue calculated from 

[S III]A18.7//m/[S III]A33.4^m. Black con- 
tours: log Ue = 2.4, 2.5, 2.6, 2.7. White contour: 
a single level of 3cm radio emission (see Figure 
The ridge is prominent, as well as a region of 
increased density to the south. 



21 



100 




-3.5 -3.0 



-1.0 -0.5 



Fig. 20. — Fits to U and Trad at a representa- 
tive location 05'*38"37" -69°06'12". Loci of con- 
sistency with the data are shown for fitting Ar, Ne, 
and S ratios alone (lines from upper left to lower 
right), the difference between the S and Ne ratios 
(log([S IV] /[S III])-log([Ne III]/[Ne II]), dashed 
lines, l-sigma confidence interval shown), and the 
combination of all ratios (bold diamond) . See text 
ijS.II for discussion. 
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Fig. 21. — Top: fitted Trad based on the S and 
Ne line ratios. Middle: fitted U. Higher values 
are dark in both grayscale images. Bottom: fit- 
ted Trad (red) and U (green), plus optical B band 
(blue) . All images are in log space except for the B 
band map. White contours: a single level of 3cm 
radio emission (see Figure[5]). The black lines indi- 
cate a scale of 1' or 15 pc. Sources of interest (see 
Tabled]) are labeled on the Trad map. The derived 
parameters. Trad and U, should be compared with 
the ratio maps in Figure [TSl 




Fig. 22. — (top) Radial dependence of conditions 
in the ionized e function of radial distance 

from R136. The different lines show azimuthal 
median values of the ionic line ratios and of the 
fitted Trad and U from photoionization models, 
(bottom) Radial parameters in the eastern and 
western directions. R136 is an off-center blister 
H II region, much closer to the bubble wall in the 
eastern than western direction. 
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Fig. 23. — Ratios of line strengths calculated from plane-parallel shock models (the actual models that were 
run are diamonds, connected for clarity). Although the Ne, Ar, and S ratios can be produced by shocks of 
a few 100km s~^ in diffuse gas, those shocks should also produce [O IV] emission, which we do not detect in 
30 Doradus. 
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6. Conclusions 

We present an infrared data cube of 30 Doradus, 
observed with the InfraRed Spectrograph and 
Multiband Imaging Spectrometer on the Spitzer 
Space Telescope. Aromatic dust emission features 
are of modest strength in the 30 Doradus region, 
concentrated in the arc-like two-lobed ridge coin- 
cident with CO emission and partially encircling 
the central R136 cluster. Detailed analysis of the 
dust content of the region from these features 
will follow in a subsequent publication, although 
the high average dust temperature is immediately 
obvious from the shape of the MIPS/SED spec- 
trum. Of the pure rotational lines of molecular 
hydrogen, only S(3)A9.67/im is detected with any 
significance, peaking on the bright "ridge" that 
dominates the morphology of the region at many 
wavelengths (see Fig [8|) . Low- ionization atomic 
lines are present but not particularly strong: [Si II] 
is detected in the outskirts of the mapped region, 
outside the region of strongest centimeter con- 
tinuum (highest emission measure ionized gas). 
[Fe II] A26.0/Ltm is only tentatively detected in two 
locations. 

Two hydrogen recombination lines Hua (7-6) 
and Hu7 (9-6) can be mapped over a large frac- 
tion of the observed area, and from these a relative 
measure of the extinction can be calculated. Inde- 
pendent estimates of the extinction are calculated 
from fitting the strength of the 10/^m silicate fea- 
ture, and from the ratio of Ha to centimeter con- 
tinuum. While none of the three extinction mea- 
sures is extremely high signal-to-noise, the loca- 
tions where all three show higher extinction should 
be quite secure - in particular three distinct loca- 
tions in the west, southwest, and northeast, out- 
side of the "ridge" (Fig [HI). 

The strongest atomic or molecular features in 
the data cube are the moderate-ionization ionic 
lines [Ar II]A7.0^m, [Ar III]A9.0/zm, [Ne II]A12.8^m, 
[Ne III]A15.5^m, [S III]A18.7/xm, [S III]A33.4^m, 
and [S IV]A10.5^m. These lines are most sensitive 
to physical conditions in the ionized gas: the ion- 
ization parameter U and hardness of the ionizing 
field, parametrized by the radiation temperature 
Trad- We fit the Neon and Sulphur ratios with 
photoionization models to derive a 2-dimensional 
map of U and Trad- We find that the excitation 
generally follows the "ridge and bubble" morphol- 



ogy, and that under the assumptions of pure pho- 
toionization, there are "hot spots" of hardened 
ionizing field in the east of the bubble and be- 
tween the two lobes of the ridge (corresponding to 
two molecular clouds). We also compare the line 
ratios to shock models and find poorer agreement 
with the data. In particular the nondetections of 
emission from more highly ionized species [O IV] 
and [Ne V] suggest that photoionization domi- 
nates over collisional excitation by shocks. Over- 
all, the local effects of hot stars in 30 Doradus 
(such as the single WR star on the eastern side 
of the bubble) appear to dominate over any large- 
scale trend with distance from the central cluster 
R136. 
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A. Flux calibration and adjustment 

A. 0.1. Background subtraction 

Our off-target baclcground observations bracket 
tfie data set in time. We chose to take an aver- 
age (witti minima and maxima trimmed) of all the 
background records and subtract that from all of 
the data records. It was possible that the two 
sets of background records could be appreciably 
different, reflecting a gradual change in the in- 
strument (e.g. electronic drift), and that a lin- 
ear interpolation of the two sets over time would 
be more appropriate as the background for each 
data record. Indeed, the mean flux level of our 
background records tends to increase with time, 
varying by 8% between the first and last expo- 
sures. Therefore, the average of the second set of 
background records has a higher mean flux level 
than the average of the first set. 

We tested three different background subtrac- 
tions on several apertures from the SLl module: 
using only the first background set, using only the 
second, and using the mean of the two sets. The 
difference between the resulting spectra was very 
small. Specifically, the difference in continuum 
levels between typical spectra using the different 
background subtractions was on the order of 2%. 
For very low signal-to-noise spectra extracted from 
faint areas of the map, the continuum difference 
could reach 10%. We emphasize that, even then, 
the shape of the spectrum and the strength of the 
emission line features did not appreciably change. 
For the highest signal-to-noise spectra, the con- 
tinuum difference was reduced to about 0.1%. We 
conclude that using the mean of all background 
records is adequate. 

A. 0.2. Flux calibration 

The IRS SL slit is approximately 3.6 arcsec- 
onds wide, and th e LL is approximately 10.6 arc- 
secon ds wide (see ISpitzer Observer's Manual 7.1 
20061 Table 7.5). Some fraction of the instru- 



ment point-spread function (PSF) falls outside 
of the slit, depending on wavelength. For ex- 
ample, the full width half max (FWHM) of the 
PSF at 14.5/L(m, at the upper wavelength range of 
the SL first order, is about 3. 7 4 arcs econds (see 



ISnitzer Observer's Manual 7.1 1 I2006L §8.1.2.1) 



flux. However, in the case of a spatially uniform 
extended source, there is no net loss of flux. CU- 
BISM can apply a slit loss correction factor which 
re-corrected the data to the original levels. In 
the case of 30 Doradus, an extended sources, we 
opted to do so. However, it is important to note 
that in some areas of the map, there are bright 
features which are neither point sources nor uni- 
formly extended emission, which have some net 
loss of flu x. The magn itude of this loss is at most 
36% Csee lShethlliooa Slide 13)0 . This loss may 
be responsible for some of the mismatches we ex- 
perienced between the flux densities in the four 
different modules at some points in the map. 

A. 0.3. Final adjustments to the spectra 

The wavelength range of the spectra needs to 
be trimmed. The full range of the IRS in each 
of the four modules includes bins where the re- 
sponse is unreliable. In the LL modules, where the 
spectrum of 30 Doradus is generally smooth, the 
drop in response at the long-wave edge of the LL2 
module and the increase in noise at the long-wave 
edge of the LLl module are evident. We trimmed 
the LL modules based on these observations. The 
spectrum of 30 Doradus is more complex in the SL 
modules, however, with many emission lines and 
an unclear continuum. In order to determine the 
range of wavelengths which can be trusted in each 
SL module, we examined the IRS low resolution 
staring mode spectrum of a Lacertae, an Al dwarf 
star whose smooth spectrum is well known. In 
each of the SL modules, the spectrum departs from 
its smooth curve at the margins. We trimmed 
our SL spectra to the wavelengths where the IRS 
spectrum of a Lac stayed faithful to the smooth 
curve. We also found, based on some spectra from 
30 Doradus, that the limit at the long-wave edge 
of the SLl module based on the spectrum of a Lac 
was not conservative enough; consequently, we de- 
creased the trusted upper limit on wavelength for 
the SLl. We could not use the spectrum of a 
Lac to determine the reliable wavelengths for the 
LL modules because the stellar spectrum has low 
signal-to-noise at such long wavelengths. 

After extracting our spectra, we found that 
some of the SL spectra were suspiciously low, even 
negative for some counts. This may be a problem 



Optimized for point sources, the Spitzer Space 
Center pipeline adjusts the data for this loss of 



* Available at: ssc.spitzer.caltech.edu/sust / workshop/2006data2 / talks /kartik.pdf 
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with background subtraction, and thus requires an 
additive correction. We raised both SL modules by 
the additive quantity necessary to make the floor 
of the SL2 spectrum (the median of the five lowest 
counts) non-negative. Over 72% of the map, this 
correction was unnecessary. Over the remaining 
28%, mainly in areas of low signal, the SL spec- 
trum was raised by a mean level of 5.5 MJy/sr, 
with standard deviation of 6.5 MJy/sr. For com- 
parison, the mean level of the background spec- 
trum in the SL is 3.2 MJy/sr and the mean level 
of the background-subtracted SL spectrum aver- 
aged across a large region of 30 Doradus is about 
100 MJy/sr. 

The different modules of an IRS spectrum do 
not always match up well in flux density. The mis- 
matches between the SLl and SL2 modules, and 
between the LLl and LL2 modules, were rarely 
serious. However, we noticed more consistent mis- 
matches between the SLl and LL2 modules. The 
primary cause is probably that the emission in the 
30 Doradus region is often best characterized as 
neither a point source nor as a uniform extended 
source, but instead as something in between. We 
have treated it as a uniform extended source (see 
i)A.0.2[) . likely losing some flux. This would call for 
a multiplicative correction. We found logarithmic 
fits to the continuum emission close to the junction 
in each module and multiplied the SL spectra by 
the ratio of the fits in the overlapping wavelengths. 
The mean value of the multiplicative correction is 
1.2 with a standard deviation of 1.2. 
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B. Quality Assurance 

In order to check our work on reducing the spec- 
tral map, we compared our spectra, taken with 
the low resolution modules of the IRS in mapping 
mode, to spectra of several areas of the 30 Doradus 
nebula taken with both the low and high resolution 
modules of the I RS, in staring mode (pri vate cor- 
respondence and Lebouteiller et al. 2007 E These 
spectra will henceforth be referred to as the GTO 
spectra. The apertures they span are displayed 
in green in Figure [3l By selecting the subsets of 
our map that coincide with the GTO apertures, 
we can produce spectra which should be consis- 



Fig. 24. — In blue, the fit to our low resolution 
spectra. In red, the same type of fit applied to the 
high resolution GTO spectra, with high resolution 
data at long wavelengths combined with low reso- 
lution data at short wavelengths (the transition is 
at about 10/.tm). Top: a relatively faint area in the 
southern part of the 30 Doradus nebula. Bottom: 
a brighter area, where the low resolution spectrum 
suffers from partial saturation and falling response 
in the LLl module. 



Sirs GTO program, ID#63, AOR keys 4382720, 12081152, 
and 12081408 
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Fig. 25. — Two comparisons of our low-resolution 
IRS data to hi K h-reso lution IRS data from 
Lebouteiller et al.l ( 2007 ). (The two panels are 



drawn from different regions on the sky, but for 
each panel the high and low resolution spectra 
are both extracted from the same aperture on 
the sky.) Both the measured spectra (symbols) 
and our fits made with PAHFIT (solid and dotted 
lines) are displayed. Note that while the high reso- 
lution spectrum has higher resolution for the lines, 
and the continuum strength may vary, the overall 
strength of the lines agree. Also note that the 
PAH has a fairly smooth appearance in the low- 
resolution plot and a more complex morphology 
in the high-resolution spectrum, which can cause 
problems when attempting to fit PAH strength in 
high-resolution. See Appendix IB] 
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Fig. 26. — Comparison of the fitted integrated 
strengths of four emission lines. On the horizontal 
axis are the results from our low resolution spec- 
tra. On the vertical axis are the results of using 
the same type of fit on the GTO spectra, which 
are high resolution longward of about 10 /im. Each 
symbol represents one area of the 30 Doradus neb- 
ula, with horizontal and vertical error bars. The 
discrepancies from the line of unity are discussed 
in Appendix |B] Among those discrepancies, the 
fitting algorithm we used is not designed to fit the 
more complex structure of PAHs that is resolved 
in the high resolution spectrum. Another reason 
for the discrepancy in the Humphreys-a is that 
our fits to the line also include the fainter H2 S(2) 
A12.28/im emission line, which is not well resolved 
at low resolution. The square symbols in the lower 
right plot call attention to two outlying points dis- 
cussed in Appendix |B] 
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tent with theirs. We appHed PAHFIT, the same 
fitting algorithm that we use on our spectral map, 
to the GTO spectra. We made some adjustments 
to PAHFIT to allow for the much higher spectral 
resolution in the GTO spectra redward of about 
10/xm, where there is coverage from the high res- 
olution modules. 

In Figure [211 the fitted spectra from two of 
these apertures are shown. The overall agree- 
ment in the shape of the spectrum is clear, with 
a notable exception in the continuum from 20 
to 32/im(likely the result of nonlinearity in the 
short- low detector response). In the upper panel, 
taken from an aperture toward the south of 
the 30 Doradus nebula, at around 05''38™50'' - 
69°06'41", where the continuum emission at long 
wavelengths is relatively faint, the agreement be- 
tween the GTO spectrum and ours is close, but the 
effect starts to become noticeable at higher flux 
densities. In the lower panel, taken from a brighter 
part of the 30 Doradus nebula at 05''38™48^ - 
69°04'10" where the falling response in the LLl 
is more severe, the difference between the GTO 
spectrum and ours is larger. 

We are most concerned with whether we can 
obtain the same emission line strengths by fitting 
the GTO spectra in the same manner as ours. Fig- 
ures [24] and [25| show that the PAHs have similar 
profiles, while the unresolved atomic lines are nar- 
rower with higher central intensities in the long- 
wave, high resolution part of the GTO spectra. 
The integrated strengths of several of the long- 
wave emission lines of interest are plotted in Fig- 
ure [26l for twelve GTO spectra and the corre- 
sponding apertures in our map. Perfect agreement 
falls along the line of unity. 

The fits to the PAH complex at 11.3/xm agree 
well only at low strength. This may be because 
of detailed fine structure that appears in the PAH 
complex at high resolution. In the low resolution 
spectrum, only a broad feature at 11.3/im, with a 
shoulder at 11.0/zm, is apparent. A Drude profile 
is sufficient to fit that, and a more complex fitting 
profile would not be appropriate because of the 
lack of resolution. In high resolution such as in 
the GTO spectra, the shoulder at 11.0/Ltm resolves 
into a completely separate feature (see Figure 1^. 
possibly another PAH feature. This might also 
be the Paschen-7 HI recombination line, although 
so significantly shifted from its rest wavelength at 



10.95/im as to make such an identification ques- 
tionable. The main PAH feature at 11.3/j,m dis- 
plays a different profile that would require a more 
complex fitting algorithm. Thus, we are not con- 
cerned about this disparity. The conclusion is that 
PAHFIT is adequate for fitting PAHs at low reso- 
lution, but may not sufficiently describe the same 
features at high resolution - not surprising, as the 
list of features used in PAHFIT was empirically 
tuned using low resolution spectra. 

A systematic disparity is evident in the com- 
parison of the Humphreys-a A12.37/im fitted line 
strength. Upon closer examination, two factors 
surfaced. Both support the trend that the ffited 
line strength at low resolution is higher than the 
line strength at high resolution. First, the fainter 
molecular hydrogen line, H2 S(2) A12.28/ini, is 
resolved in high resolution. The fit to the high 
resolution Humphreys-a line therefore does not 
include any of the flux from the molecular hy- 
drogen line. Meanwhile, in the low resolution 
spectrum, the flux from both lines contributes to 
the same broad feature. We are confident that 
Humphreys-a contributes most of the flux wc see 
in low resolution, in part because the line cen- 
ter more closely coincides with 12.37 than with 
12.28/im, and in part because, as we can see in 
the high resolution spectra, the Humphreys-a line 
is equal to or stronger than H2 S(2) in all of the 
areas of our map covered by the high resolution 
spectra. However, the flux of both emission lines 
is contributing to the low-resolution ffi, and so to 
some extent, our fits to the line strength of the 
Humphreys-a line are overestimated. Quantita- 
tive comparison between high- and low-resolution 
spectra and between fits with and without the 
molecular hydrogen line show that the contamina- 
tion of the Humphreys-a line strength is at most 
~20% (see §4.ip . This effect is weakest in the parts 
of 30 Doradus dominated by ionized gas that are 
the focus of this paper, but should be kept in mind 
for analysis of the outer parts of the region and 
PDR physics in future papers. The second fac- 
tor that contributes to the disparity seen in the 
second panel of Figure [26] is the PAH feature at 
12.6/im. As discussed above in the context of the 
PAH complex at 11.3/im, the fitting package PAH- 
FIT is intended for the relatively smooth features 
seen in low resolution, and does not handle the 
complex morphology of PAHs seen at high resolu- 
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tion as well. In the case of the PAH at 12.6/xm, 
its line strength in the high resolution fit is gener- 
ally overestimated. Its broad wings thus lead the 
strength of the Humphreys-a line to be underesti- 
mated in the high resolution fit. The slight overes- 
timation of the line in the low resolution spectrum 
and the slight underestimation of the line in the 
high resolution spectrum both contribute to the 
disparity. 

The results for the forbidden neon lines, as for 
most of the stronger lines, show a tighter corre- 
lation. The close agreement in the fits to the 
low and high resolution spectra for the weaker 
[Ne II]A12.8/im is particularly encouraging be- 
cause one of our concerns in interpreting the neon 
ratio was the possible entanglement of the [Ne II] 
atomic emission line with the above-mentioned 
PAH at 12.6/im. In the high resolution spectrum, 
the PAH (with its often-unusual shape) and the 
atomic line are resolved. As in the case of the 
Humphreys-a line, the strength of the [Ne II] line 
may be slightly underestimated because of the 
poor, overestimated fit to the PAH at 12.6/im. 
The close agreement between the low resolution 
and high resolution fit, across almost all of the 
twelve apertures, regardless of PAH strength in 
each aperture, indicates that there are no system- 
atic errors to the low resolution fit that do not also 
appear in the high resolution fit. Thus, our results 
may share the slight bias toward underestimating 
the [Ne II] strength, but we are generally able to 
decompose [Ne II] and the PAH at 12.6//m despite 
the lack of spectral resolution. 

Finally, the correlation for the fitted line 
strength of [Ne III]AI5.5/im is similarly good, with 
the exception of two apertures where the high res- 
olution fitted line strength is significantly higher 
than the low resolution fitted line strength. These 
two apertures (displayed with squares rather than 
diamonds in Figure [26]) are located in a small re- 
gion of the map that exemplifies the extreme of 
the trend shown in Figure [24l where the response 
in the LLI module falls until the module finally 
saturates. In these cases, our fit to the continuum 
with thermal dust components is not very success- 
ful. The fitted continuum is dragged down in the 
area from 20 to 32/im, which forces it up in the 
region of I5/im. This causes our estimation of the 
strength of [Ne III] to be underestimated. 

For these two anomalous apertures, we con- 



ducted a test by removing the faulty LLI spec- 
trum and fitting only the remaining spectrum, to 
get a better estimation of continuum in the re- 
gion of [Ne III]. This increased the fitted strength 
of [Ne III], though not enough to get it to agree 
with the value we obtained by fitting the high- 
resolution spectrum. 

The effect of long- wave line strengths being un- 
derestimated in our low resolution spectra com- 
pared to the high resolution GTO spectra is seen 
for [Ne III]AI5.5/im and [S III]A33.4/im in some 
areas of the map. Empirically, the line strengths 
of [Ne II]A12.8/im and [S III]AI8.7/im appear to 
be unaffected, as shown by their close agreement 
with the GTO spectra. It is crucial to note 
that this effect is restricted to the small areas 
of the map where we can see falling response 
in the LLI module, as drawn in Figure SI In 
that localized area, we can expect our values of 
[Ne III] and [S III]A33.4/im to be underestimated 
relative to the GTO results by perhaps 30%, as 
shown in Figure [26] (the effect is very similar for 
[S III]A33.4/im). Outside these small areas of the 
map, the agreement between our spectra and the 
high-resolution GTO spectra is excellent for these 
two emission lines. 
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